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Abstract — We  consider  ad  hoc  wireless  networks  that  use  directional  antennas  and  have  limited  energy  resources.  To  explore 
quantitatively  the  advantage  offered  by  the  use  of  directional  antennas  over  the  case  of  omnidirectional  antennas,  we  consider  the 
case  of  connection-oriented  multicast  traffic.  Building  upon  our  prior  work  on  multicasting  algorithms,  we  introduce  two  protocols  that 
exploit  the  use  of  directional  antennas  and  evaluate  their  performance.  We  observe  significant  improvement  with  respect  to  the 
omnidirectional  case,  in  terms  of  both  energy  efficiency  and  network  lifetime.  Additionally,  we  show  that  further  substantial  increase  in 
the  network’s  lifetime  can  be  achieved  by  incorporating  a  simple  measure  of  a  node's  residual  energy  into  the  node's  cost  function. 

Index  Terms — Broadcast,  multicast,  energy  efficient,  directional  antenna,  ad  hoc  network. 
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1  Introduction 

he  use  of  directional  antennas  can  provide  energy  savings 
and  interference  reduction  by  concentrating  RF  energy 
where  it  is  needed.  Hence,  they  are  especially  useful  in 
networks  with  finite  energy  resources  and  a  limited  number 
of  available  frequencies.  In  this  paper,  we  develop  and 
evaluate  algorithms  for  broadcasting  and  multicasting  that 
are  suitable  for  use  in  networks  with  directional  antennas  and 
limited  battery  capability,  and  compare  performance  to  that 
achieved  when  antennas  are  omnidirectional.  We  focus  on 
the  problem  of  tree  construction  for  source-initiated,  session- 
based  traffic  in  all-wireless  (i.e.,  infrastructureless,  peer-to- 
peer,  or  ad  hoc)  multihop  networks. 

In  our  earlier  studies  [1],  [2],  we  developed  energy-aware 
algorithms  for  the  construction  of  broadcast  and  multicast 
trees  for  networks  with  omnidirectional  antennas.  In  the 
broadcasting  case,  the  goal  is  to  find  a  minimum-cost  (e.g., 
minimum  power)  tree,  rooted  at  the  Source  node,  that 
reaches  all  other  nodes  in  the  network;  in  the  multicasting 
case,  a  specified  subset  of  the  nodes  must  be  reached,  while 
others  may  serve  as  relays.  These  algorithms  are  known  as 
Broadcast  Incremental  Power  (BIP)  and  Multicast  Incre¬ 
mental  Power  (MIP)  [1],  [2],  respectively.1  It  has  subse¬ 
quently  been  proven  that  the  minimum-cost  broadcasting 

1.  BIP  has  also  been  studied  by  other  researchers  [3],  [4],  and  [5]. 
Alternative  algorithms,  which  perform  better  than  BIP/MIP  but  at  the 
expense  of  higher  complexity,  were  introduced  in  [6]. 
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problem  in  wireless  networks  with  omnidirectional  anten¬ 
nas  is  NP-complete  [7],  [4].  By  contrast,  the  minimum-cost 
broadcasting  problem  in  wired  networks  can  be  formulated 
as  the  well-known,  and  easily  solved,  minimum-cost 
spanning  tree  (MST)  problem. 

BIP  and  MIP  are  examples  of  "node-based"  algorithms, 
which  exploit  the  "wireless  multicast  advantage"  property 
associated  with  omnidirectional  antennas,  namely,  the 
capability  for  a  node  to  reach  several  neighbors  by  using 
a  transmission  power  level  sufficient  to  reach  the  most 
distant  one.  This  property  no  longer  holds  for  directional 
antennas  because  the  RF  power  needed  to  sustain  sufficient 
signal  levels  depends  on  the  beamwidth.  In  this  paper, 
using  the  incremental  power  philosophy  as  a  starting  point, 
we  demonstrate  the  issues  that  arise  when  directional 
antennas  are  used,  we  develop  algorithms  for  tree  con¬ 
struction,  and  we  illustrate  trade-offs  between  complexity 
and  performance. 

We  address  energy-aware  networking  from  the  per¬ 
spectives  of  both  energy-efficient  and  energy-limited  opera¬ 
tion.  In  the  former  case,  it  is  implicitly  assumed  that  energy 
reserves  (batteries)  can  be  renewed  during  the  course  of 
network  operation,  and  a  cost  is  associated  with  the 
quantity  of  energy  that  is  expended.  In  the  latter  case,  each 
node  is  equipped  with  batteries  that  cannot  be  recharged 
during  network  operation.  Thus,  there  is  a  hard  constraint 
of  a  fixed  quantity  of  energy  at  each  of  the  network  nodes.  We 
briefly  address  some  of  the  fundamental  differences 
between  energy-limited  and  energy-efficient  network  op¬ 
eration.  Our  primary  focus  in  this  paper  is  on  energy- 
limited  operation. 

In  [8],  for  the  case  of  omnidirectional  antennas,  we 
studied  the  performance  of  MIP  under  the  constraint  of  a 
fixed  quantity  of  energy  at  each  of  the  network  nodes.  We 
demonstrated  that  the  lifetime  of  the  network  can  be 
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extended  significantly  by  incorporating  into  the  tree- 
construction  process  a  cost-function  that  reflects  the 
residual  energy  at  the  nodes.  Here,  we  demonstrate  that, 
when  modified  to  exploit  the  properties  of  directional 
antennas,  the  incremental  power  philosophy  of  BIP  and 
MIP  can  provide  further  substantial  gain  in  both  energy 
efficiency  and  network  lifetime. 

Few  studies  of  the  use  of  directional  antennas  for  energy- 
aware  operation  have  appeared  in  the  literature.  Our 
preliminary  results  on  the  work  presented  in  this  paper 
appeared  in  [9],  and  the  complete  model  was  presented  in 
[10].  The  use  of  directional  antennas  to  provide  improved 
energy  efficiency  and  interference  reduction  in  multicast 
applications  (but  assuming  that  the  trees  are  already 
known)  was  later  studied  in  [11].  Other  applications  of 
the  use  of  directional  antennas  for  energy  efficiency  were 
studied  in  [12]. 

The  focus  of  this  paper  is  the  construction  of  energy- 
aware  broadcast  or  multicast  trees,  specifically  for  the  case 
of  directional  antennas.  To  provide  a  better  understanding 
of  the  fundamental  issues  involved,  as  in  [1],  [2],  [8],  we  do 
not  address  many  of  the  implementation  issues  that  will 
have  to  be  addressed  eventually.  For  example,  we  do  not 
address  admission  control;  we  assume  that  a  tree  is 
constructed  to  reach  all  destinations  that  are,  in  fact, 
reachable,  regardless  of  the  cost  (energy)  to  do  so.  Nor  do 
we  address  the  control  of  transmission  rate  or  other 
transmission  parameter  (which  can  affect  session  duration, 
energy  usage,  quality  of  service,  etc.).  We  also  choose  to 
assume  that  the  channel  bandwidth  and  signal  design 
parameters  are  set  so  that  the  bit  rate  is  fixed.  Additionally, 
we  do  not  address  media-access  issues  or  mobility.  Thus, 
our  study  of  tree  construction  in  the  wireless  environment 
translates  to  choosing  transmission  power  and  the  set  of 
receiving  neighbor  nodes  at  each  level  in  the  multicast  tree. 

In  Section  2,  we  provide  a  brief  discussion  of  the 
differences  between  energy-constrained  and  energy-effi¬ 
cient  communication.  In  Section  3,  we  discuss  the  commu¬ 
nication  model,  including  the  impact  of  the  use  of 
directional  antennas  on  energy  consumption.  In  Section  4, 
we  define  the  multicasting  problem,  including  the  perfor¬ 
mance  measures  used  in  this  study.  In  Section  5,  we  present 
our  algorithms  for  broadcast  and  multicast-tree  construc¬ 
tion,  which  exploit  the  properties  of  directional  antennas.  In 
Section  6,  we  show  how  the  impact  of  resource  limitations 
(transceivers,  frequencies,  and  energy)  can  be  incorporated 
into  our  model.  In  Section  7,  we  present  our  performance 
results  for  energy-constrained  systems  and  demonstrate 
that  significant  increase  in  delivered  traffic  volume  can  be 
achieved  by  using  algorithms  that  exploit  the  properties  of 
directional  antennas.  Finally,  in  Section  8,  we  present  our 
conclusions  from  this  research. 

2  Energy-Constrained  vs.  Energy-Efficient 
Communication 

The  introduction  of  hard  constraints  on  the  total  amount  of 
energy  available  at  each  node  results  in  a  problem  that  is 
very  different  from  that  in  which  unlimited  energy  is 
available  (although  energy  efficiency  still  may  be  desired). 


Under  such  hard  constraints  on  energy,  the  network  is 
capable  of  operation  for  a  limited  period  of  time.  A  node 
dies  (and,  hence,  can  no  longer  transmit)  when  its  energy  is 
depleted,  and  the  network  dies  when  it  is  no  longer  capable 
of  providing  a  minimum  acceptable  level  of  service.  By 
contrast,  when  the  goal  is  energy  efficiency  (e.g.,  delivering 
the  largest  number  of  bits  per  unit  energy),  it  is  implicitly 
assumed  that  ample  energy  is  available;  in  such  cases,  the 
use  of  energy  is  essentially  treated  as  a  cost  function. 

Energy-efficient  operation  does  not  ensure  good  perfor¬ 
mance  in  energy-constrained  applications.  For  example,  use 
of  the  most  energy-efficient  routes  (or  multicast  trees)  may 
result  in  premature  depletion  of  energy  at  some  nodes.  We 
use  the  term  energy-aware  to  refer  to  systems  designed  with 
either  of  these  criteria  in  mind. 

A  more  complete  discussion  of  energy-limited  versus 
energy-efficient  communication  is  provided  in  [8],  where 
performance  results  are  provided  for  networks  that  use 
omnidirectional  antennas. 

3  The  Model 

We  consider  source-initiated,  circuit-switched,  multicast 
sessions.  The  maintenance  of  a  session  requires  the 
dedication  of  a  transceiver  at  each  participating  node 
(source  node,  relay  nodes,  and  destination  nodes),  as  well 
as  the  needed  amount  of  interference-free  bandwidth, 
throughout  the  duration  of  the  session.  The  network 
consists  of  N  nodes,  which  are  randomly  distributed  over 
a  specified  region.  Each  node  has  T  transceivers  and  can 
thus  support  up  to  T  multicast  sessions  simultaneously.  We 
assume  that  there  is  a  total  of  F  frequencies  available  to  the 
network.  Frequencies  can  be  reused,  provided  that  doing  so 
does  not  create  interference.  Thus,  congestion  (and,  hence, 
the  inability  to  reach  one  or  more  destinations)  may  arise 
when  either  an  insufficient  number  of  transceivers  or  an 
insufficient  number  of  frequencies  are  available  at  one  or 
more  nodes  along  the  path.  Alternatively,  energy-inefficient 
paths  may  have  to  be  used  when  the  best  paths  are  not 
available. 

It  is  also  of  interest  to  study  systems  that  use  time- 
division  multiple  access  (TDMA),  rather  than  multiple 
transceivers  or  multiple  channels,  to  support  multiple 
sessions  simultaneously.  In  TDMA-based  systems,  the  need 
to  assign  specific  time  slots  creates  a  much  more  difficult 
problem  than  that  of  simply  assigning  any  transceiver  (of 
perhaps  several  available)  to  a  new  session.  Alternatively,  it 
would  be  possible  to  consider  code-division  multiple  access 
(CDMA)  [13].  The  study  of  TDMA-  and  CDMA-based 
systems  is  not  pursued  here  since  we  want  to  place 
emphasis  on  the  energy  constraint  with  as  little  complica¬ 
tion  from  the  MAC  layer  as  possible. 

Any  node  is  permitted  to  initiate  multicast  sessions. 
Multicast  requests  and  session  durations  are  generated 
randomly  at  the  network  nodes.  Each  multicast  group 
consists  of  the  source  node  plus  at  least  one  destination 
node.  Additional  nodes  may  be  used  as  relays  either  to 
provide  connectivity  to  all  members  of  the  multicast  group 
or  to  reduce  overall  energy  consumption.  The  set  of  nodes 
that  support  a  multicast  session  (the  source  node,  all 
destination  nodes,  and  all  relay  nodes)  is  referred  to  as  a 
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Fig.  1 .  Use  of  directional  antenna  to  reach  two  neighboring  nodes. 

multicast  tree.  Notice  the  difference  between  this  definition 
and  the  conventional  one  that  is  based  on  links  (or  edges); 
here,  the  links  are  incidental  and  their  existence  depends  on 
the  transmission  power  of  each  node.  Thus,  it  is  the  nodes 
(rather  than  the  links)  that  are  the  fundamental  units  in 
constructing  the  tree. 

The  connectivity  of  the  network  depends  on  the 
transmission  power  and  antenna  pattern.  We  assume  that 
each  node  can  choose  its  RF  power  level  pRF,  such  that 
Pmin  <  PRF  <  Pmax-  The  nodes  in  any  particular  multicast 
tree  do  not  necessarily  have  to  use  the  same  power  levels; 
moreover,  a  node  may  use  different  power  levels  for  the 
various  multicast  trees  in  which  it  participates. 

3.1  Propagation  Model 

When  considering  omnidirectional  antennas  and  uniform 
propagation  conditions,  we  assume  that  the  received  signal 
power  is  equal  to  pr~a,  where  p  is  the  transmission  power,  r 
is  the  distance,  and  a  is  a  parameter  that  typically  takes  on  a 
value  between  two  and  four,  depending  on  the  character¬ 
istics  of  the  communication  medium.  Based  on  this  model, 
the  transmitted  power  required  to  support  a  link  between 
two  nodes  separated  by  distance  r  is  proportional  to  r"  since 
the  received  power  must  exceed  some  threshold.2  Let  us 
first  consider  the  case  of  point-to-point  communication 
between  Node  i  and  Node  j.  Without  loss  of  generality,  we 
set  the  proportionality  constant  equal  to  1,  resulting  in: 

pRF  =  RF  power  needed  for  link  between  Nodes  i  and  j 

13  (1) 
—  max{rl;.  Pmin\  ■ 

where  ry  is  the  distance  between  Node  i  and  Node  j.  The 
use  of  a  nonzero  value  of  pmm  is  a  way  to  account  for  the  fact 
that  the  r~a  dependence  applies  only  in  the  far-field  region 
(i.e.,  even  when  two  nodes  are  arbitrarily  close  to  each 
other,  a  nonzero  power  level  pm jn  is  required  to  support 
communication  between  them).  As  a  consequence  of  the 
“wireless  multicast  advantage"  property  of  omnidirectional 
systems  [14],  all  nodes  whose  distance  from  Node  i  does  not 
exceed  rtJ  will  be  able  to  receive  the  transmission  with  no 
further  energy  expenditure  at  Node  i. 

2.  This  threshold  depends  on  factors  such  as  signal  parameters,  detector 
structure,  and  noise  levels  (including  other-user  interference).  In  this  paper, 
we  assume  that  these  characteristics  are  fixed;  thus,  the  required  level  of 
received  power  is  the  same  at  all  nodes.  Thus,  we  neglect  fading  effects  that 
arise  in  wireless  channels. 


The  use  of  directional  antennas  can  permit  energy  savings 
and  reduce  interference  by  concentrating  transmission 
energy  where  it  is  needed.  On  the  other  hand,  since  only  the 
nodes  located  within  the  transmitting  node's  antenna  beam 
can  receive  the  signal,  the  effect  of  the  wireless  multicast 
advantage  may  be  diminished;  additional  power  is  needed  as 
the  beamwidth  is  increased  to  include  additional  receiving 
nodes.  We  use  an  idealized  model  in  which  we  assume  that  all 
of  the  transmitted  energy  is  concentrated  uniformly  in  a  beam 
of  width  9,  as  shown  in  Fig.  1,  where  Node  i  transmits  to  Node 
j  and  Node  k.3  Then,  the  RF  power  needed  by  a  node  to 
transmit  to  a  distance  r  using  beamwidth  9  (and  thus  to  reach 
all  nodes  in  this  sector)  is 

pRF(r,  9)  =  max  jr“  ^  ,pmin|  •  (2) 

Consequently,  the  use  of  narrow  beams  permits  energy 
saving  (for  a  given  communication  range)  or  range 
extension  (for  a  given  transmitter  power  level),  as  compared 
to  the  use  of  omnidirectional  antennas.  Specifically,  for  a 
given  value  of  pmax,  the  maximum  range  is  increased  by  a 
factor  of  (360/f?)1^0,  compared  with  the  case  of  omnidirec¬ 
tional  antennas. 

We  assume  that  the  beamwidth  6  can  be  chosen  so  that 
#min  <  9  <  #max.  Furthermore,  we  assume  that  each  node 
knows  the  precise  locations  of  its  potential  neighbors  (i.e., 
all  nodes  that  can  be  reached  with  acceptable  signal-to- 
interference  ratio  when  transmitted  with  pRF  <  pmax)  and 
that  each  antenna  beam  can  be  pointed  in  any  desired 
direction  to  provide  connectivity  to  a  subset  of  the  nodes 
that  are  within  communication  range.  (In  practice,  the 
number  of  antenna  elements  needed  tends  to  increase  as 
9m in  decreases.) 

We  assume  that  one  antenna  beam  can  be  supported  for 
each  session  in  which  a  node  participates;  thus,  the  use  of 
directional  antennas  does  not  have  an  impact  on  the 
number  of  sessions  that  a  node  can  support  simultaneously 
(as  compared  to  an  implementation  with  omnidirectional 
antennas).  Additionally,  both  9  and  the  direction  in  which 
the  beam  points  are  chosen  independently  for  each  session 
in  which  a  node  participates.  Although  setting  9  =  0min  is 
appropriate  for  point-to-point  applications,  it  is  often 
appropriate  to  use  larger  values  of  9  in  multicast  applica¬ 
tions  since  a  node  may  have  several  downstream  neighbors, 
all  of  which  must  be  included  in  a  single  beam  (based  on  the 
assumption  just  made  above).  We  discuss  the  choice  of  9  in 
our  discussion  of  multicast  algorithms  in  Section  5. 

Although  we  do  consider  energy  expenditures  asso¬ 
ciated  with  processing  at  each  node  (in  addition  to  that  for 
RF  transmission),  we  do  not  explicitly  connect  the  amount 
of  processing  energy  with  the  beamwidth  of  the  antenna. 
This  coupling  is  deferred  for  future  investigation. 

The  use  of  directional  receiving  antennas  would  also 
have  a  beneficial  impact  since  background  noise  and  other- 
user  interference  would  be  troublesome  only  when  located 
within  the  antenna  beamwidth  rather  than  the  entire 
omnidirectional  region.  Thus,  lower  signal  levels  would 
be  needed  to  provide  the  required  performance.  However, 

3.  In  realistic  systems,  antenna  patterns  would  be  characterized  by 
nonuniform  gain  in  the  primary  sector  of  coverage,  as  well  as  by  sidelobes. 
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we  assume  the  use  of  omnidirectional  receiving  antennas  to 
simplify  the  model. 

It  is  also  possible  to  consider  alternative  models,  which 
may  incorporate  one  or  more  of  the  following: 

•  fixed  beamwidth  (i.e.,  0m;n  =  9max), 

•  a  single  beam  per  node, 

•  multiple  beams  per  session, 

•  constraint  on  number  of  beams  per  node  (possibly  >  T), 
and 

•  directional  receiving  antennas. 

However,  these  are  not  addressed  in  this  paper. 

3.2  Energy  Expenditure 

In  addition  to  RF  propagation,  energy  is  also  expended  for 
transmission  (encoding,  modulation,  etc.)  and  reception 
(demodulation,  decoding,  etc.).  We  define: 

•  pT  =  transmission  processing  power  and 

•  pR  =  reception  processing  power. 

We  assume  that  these  quantities  are  the  same  at  all 
nodes,  and  we  neglect  any  energy  consumption  occurring 
when  the  node  is  simply  "on"  without  transmitting  or 
receiving,  although  it  would  be  straightforward  to  include 
it.  However,  for  networks  with  regularly  steady  traffic,  it  is 
rare  that  a  node  will  not  be  receiving  or  transmitting  at  any 
given  time,  especially  in  view  of  our  assumption  that 
multiple  transceivers  are  available  at  each  node.  The  total 
power  expenditure  of  a  node,  when  transmitting  to  a 
maximum  range  r  over  a  sector  of  width  9,  is 

p  =  pRF{r,  6)  +  pT  +  pR  l(node  is  not  the  source),  (3) 

where  pRF(r ,  6)  is  defined  in  (2),  and  the  indicator  function 
is  included  because  the  pR  term  is  not  needed  for  the  source 
node.  A  leaf  node,  since  it  does  not  transmit  but  only 
receives,  has  a  total  power  expenditure  of  pR. 

We  assume  that  each  node  starts  with  a  finite  quantity  of 
battery  energy.4  For  example,  Node  i  has  energy  £)(0)  at 
time  0.  The  residual  energy  at  Node  i  at  time  t  is 

Ei(t)  =  £)(0)  -  T  Pi(r)dT,  (4) 

Jo 

where  Pi(r)  is  the  total  power  expended  at  Node  i  at  time 
r.5  We  say  that  a  node  is  "alive"  as  long  as  its  residual 
energy  is  positive  and  that  it  dies  when  its  residual  energy 
decreases  to  zero.  Based  on  our  assumptions,  a  "dead"  node 
cannot  participate,  even  as  a  receive-only  leaf  node. 

4  The  Multicasting  Problem 

The  establishment  of  a  multicast  tree  requires  the  specifica¬ 
tion  of  the  transmitted  power  levels,  the  frequencies  used 
by  each  node,  and  the  commitment  of  the  needed 

4.  We  assume  that  the  battery  has  a  fixed  capacity,  i.e.,  we  neglect  the  fact 
that  the  total  energy  that  can  be  supplied  by  a  battery  depends  in  part  on 
the  discharge  rate  and  duty  cycle  [15].  We  also  neglect  any  nonlinear 
behavior,  which  may  characterize  power  amplifiers  especially  at  high 
output  levels. 

5.  Since  Node  i  may  be  transmitting  as  a  member  of  several  trees 
simultaneously,  Pi(r)  is  the  sum  of  the  powers  for  all  such  trees  at  time  r. 


transceiver  resources  throughout  the  duration  of  the 
session. 

We  assume  that  multicast  session  requests  arrive  to  each 
of  the  N  nodes  at  rate  \/N  arrivals  per  unit  time.  The  set  of 
desired  destinations  is  chosen  randomly  for  each  arrival. 
We  say  that  a  destination  can  be  reached  if  the  following 
conditions  are  satisfied: 

•  there  exists  a  path  from  the  source  to  it  (i.e.,  the 
transmitted  power  required  to  support  the  path  does 
not  exceed  pmax  at  any  node), 

•  a  transceiver  is  available  (i.e.,  not  already  supporting 
another  session)  at  each  node  along  the  path,  and 

•  a  suitable  frequency  assignment  can  be  found  to 
support  the  path  (i.e.,  a  noninterfering  frequency  is 
available  to  support  the  link  between  each  node  pair 
in  the  network  along  the  path;  these  frequency 
assignments  must  not  interfere  with,  or  suffer 
interference  from,  currently  ongoing  sessions). 

As  noted  earlier,  all  multicast  requests  are  accepted  as 
long  as  one  or  more  of  the  intended  destinations  can  be 
reached  and  paths  are  established  to  all  reachable  destina¬ 
tions,  regardless  of  the  cost  required  to  do  so. 

4.1  Performance  Measures 

In  this  paper,  we  focus  primarily  on  one  particular 
performance  measure,  which  is  especially  well  suited  for 
energy-limited  applications,  namely,  the  total  delivered 
traffic  volume  during  the  lifetime  of  the  network.  We  also 
consider  the  related  quantity  of  traffic  volume  per  unit 
energy. 

We  first  introduce  some  notation.  We  assume  that,  once  a 
session  (multicast  tree)  is  established,  communication  takes 
place  at  a  constant  rate  of  R  bits/s,  which  is  the  same  for 
each  session  request,  and  which  is  independent  of  A. 
Session  duration  is  exponentially  distributed  with  mean 

1/p  =  1. 

Since  partial  multicast  sessions  may  take  place  (because 
some  nodes  may  be  unreachable),  the  performance  metric 
should  provide  a  reward  that  reflects  the  number  of 
destinations  that  are  actually  reached.  We  define 

•  Hi  =  #  of  intended  destinations  in  session  i, 

•  m,  =  #  of  destinations  reached  in  session  i, 

•  Pi  =  sum  of  the  transmitter  powers  used  by  all  nodes 
in  session  i,  and 

•  di  =  duration  of  session  i. 

4. 1. 1  Delivered  Traffic  Volume 

The  delivered  traffic  volume  is  directly  proportional  to  both 
the  number  of  destinations  that  are  reached  and  to  the 
duration  of  each  session.  Specifically,  each  destination  node 
participating  in  multicast  session  i  receives  hi  =  R  di  bits 
during  the  course  of  the  session.  The  total  quantity  of  data 
delivered  during  session  i  is  then 

B,  =  total  number  of  bits  delivered  to  all 
reached  destinations  in  session  i 
=  771,6,;. 
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Then,  the  total  quantity  of  information  delivered  to  all 
destinations  over  an  observation  interval  of  X  multicast 
requests  is: 

x  x 

Btotai  =  J2bi  =  R.J2  midi.  (5) 

i=  1  i—1 

4. 1.2  Delivered  Traffic  Volume  per  Unit  Energy 

The  energy  expenditure  in  session  i  is  Pidi.  Thus,  the  total 

energy  expenditure  over  the  observation  interval  is 

Ex  =  YJPrdl-  (6) 

i—l 

Therefore,  the  delivered  traffic  volume  per  unit  energy  over 
an  interval  of  X  arrivals  is 

r>  _BT  _RJ2hmldl 


4.1.3  Local  Cost  Metrics 

Tree  formation  consists  of  the  specification  of  transmitting 
nodes  and  their  downstream  neighbors.  When  omnidirec¬ 
tional  antennas  are  used,  it  is  sufficient  to  specify  the  set  of 
transmitting  nodes  and  their  RF  transmission  power  levels; 
when  directional  antennas  are  used,  the  antenna  pattern 
must  also  be  specified.  It  is  not  feasible  to  find  the  multicast 
trees  that  guarantee  the  optimal  values  of  global  perfor¬ 
mance  measures  such  as  Btotal,  etc.  Therefore,  we  have 
focused  on  the  development  of  "local"  strategies  that 
depend  on  "local"6  metrics,  which  find  the  multicast  tree 
that  attempts  to  minimize  an  appropriate  cost  function  for 
each  new  multicast  request. 

In  particular,  the  basic  approach  taken  in  [1]  and  [2]  is  to 
minimize  the  power  needed  to  maintain  the  tree  associated 
with  each  newly  arriving  session.7  This  power  includes  the 
RF  transmission  power  of  all  transmitting  nodes  as  well  as 
the  signal  processing  power  expended  at  transmitting  and 
receiving  nodes.  We  recognize  that  local  optimization  does 
not  guarantee  global  optimization,  e.g.,  minimizing  tree 
power  does  not  guarantee  the  minimization  of  energy  over 
an  observation  interval  of  many  arrivals.  Moreover,  even  if 
it  were  possible  to  do  so,  this  would  certainly  not  guarantee 
the  optimization  of  the  desired  global  performance  mea¬ 
sures.  Nevertheless,  it  has  been  our  experience  that  this 
approach  works  reasonably  well. 

The  multicasting  problem  is  similar  to  the  broadcasting 
problem,  except  that  only  a  specific  subset  of  the  nodes  are 
required  to  be  in  the  tree.  It  is  well-known  that  the 
determination  of  a  minimum-cost  multicast  tree  in  wired 
networks  is  a  difficult  problem,  which  can  be  modeled  as 
the  NP-complete  Steiner  tree  problem,  even  though  the 
broadcasting  problem  is  easily  formulated  as  the  MST 
problem,  which  has  low  complexity.  In  wireless  networks, 
even  the  broadcasting  problem  is  NP-complete  [7],  [4]. 
Thus,  heuristics  are  needed  for  both  broadcasting  and 

6.  "Global"  is  used  here  to  refer  to  optimization  over  a  long  observation 
interval.  "Local"  is  used  here  both  in  the  sense  of  time-local  (i.e.,  for  each 
arrival  of  a  multicast  session  request),  as  well  as  in  the  topological  sense 
(i.e.,  pertaining  to  an  individual  link  or  node). 

7.  In  Section  6,  we  introduce  a  cost  metric  that  also  involves  the  residual 

energy  at  each  node. 


multicasting.  The  two  basic  approaches  we  have  used  for 
multicasting  are  the  "pruning"  of  broadcast  trees  and  the 
superposition  of  unicast  paths  [1],  [2]. 

5  Algorithms  for  Broadcasting  and 
Multicasting  with  Directional  Antennas 

We  have  considered  two  basic  approaches  for  broadcasting 
and  multicasting  with  directional  antennas: 

•  Construct  the  tree  by  using  an  algorithm  designed 
for  omnidirectional  antennas;  then,  reduce  each 
antenna  beam  to  the  minimum  possible  width  that 
can  support  the  tree. 

•  Incorporate  directional  antenna  properties  into  the 
tree-construction  process. 

The  first  approach  can  be  used  with  any  tree-construction 
algorithm.  The  "beam-reduction"  phase  is  performed  after 
the  tree  is  constructed  by  using  an  additional  "postproces¬ 
sing"  algorithm,  which  is  appended  to  the  tree-construction 
algorithm.  The  second  approach,  which  requires  decisions 
on  beamwidth  to  be  made  at  each  step  of  the  tree 
construction  process,  can  be  used  only  with  algorithms 
that  construct  trees  by  adding  one  node  at  a  time,  such  as 
BIP  (and  its  multicasting  counterpart  MIP).  In  this  section, 
we  describe  these  approaches  in  detail. 

5.1  An  Approach  Based  on  Beamwidth-Reduction: 
Reduced  Beam  BIP  (RB-BIP)  and  MIP  (RB-MIP) 

First,  a  low-cost  broadcast  or  multicast  tree  is  formed,  using 
any  tree-construction  algorithm  (e.g.,  BIP  or  MIP),  under  the 
assumption  that  the  transmitting  antennas  are  omnidirec¬ 
tional.  Then,  after  the  tree  is  constructed  in  this  manner, 
each  transmitting  node's  antenna  beamwidth  is  reduced  to 
the  smallest  possible  value  that  provides  coverage  of  the 
node's  downstream  neighbors,  subject  to  the  constraint 
0mm  <  9  <  360.  Thus,  the  tree  structure  is  independent  of 
0mm-  We  assume  perfect  antenna  patterns  that  provide 
uniform  gain  throughout  the  cone  of  beamwidth  6  (with  no 
sidelobes),  so  it  is  not  necessary  to  extend  6  beyond  the 
direction  of  the  nodes  at  the  edges  of  the  cone.  When 
applied  to  BIP,  the  resulting  scheme  is  called  Reduced- 
Beam  BIP  (RB-BIP);  when  applied  to  MIP,  the  resulting 
scheme  is  called  RB-MIP. 

As  noted  in  Section  3.1,  we  assume  that  each  node  knows 
the  location  of  all  of  its  potential  neighbors  (at  pRF  <  Pmax), 
and  can  thus  point  its  beam  in  the  appropriate  direction  to 
provide  connectivity  to  all  of  its  children  in  the  tree.  Fig.  2a 
shows  Node  s  with  its  children  (as  determined  by  BIP). 
Once  the  set  of  children  is  determined,  the  communication 
range  r  is  the  distance  to  the  most-distant  child.  To 
minimize  the  required  power  (for  a  given  range),  the 
minimum  beamwidth  that  reaches  all  children  must  be 
determined.8  This  problem  involves  finding  the  largest 
sector  that  can  be  omitted  from  omnidirectional  coverage, 
while  still  reaching  all  of  the  children.  In  Fig.  2a,  the  choice 
of  the  appropriate  sector  of  coverage  (which  excludes  the 
sector  of  width  360  —  9)  is  obvious  from  the  graphical 

8.  In  some  applications,  it  may  be  desirable  to  incorporate  interference 
considerations  into  the  choice  of  the  sector  that  is  to  be  omitted;  however, 
we  do  not  do  so  in  this  paper. 
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(b) 


Fig.  2.  Choosing  the  sector  to  be  covered  by  the  beam,  (a)  Sector  of  beamwidth  9  and  (b)  angles  used  to  determine  sector. 


representation.  However,  it  is  helpful  to  discuss  the 
algorithmic  procedure  used  to  determine  this  sector. 

Given  the  coordinates  of  the  nodes,  we  would  like  to 
determine  the  sector  of  minimum  beamwidth  that  includes 
all  children  nodes.  In  general.  Node  s  has  n  children,  which 
are  a  subset  of  the  nodes  in  {1, 2, . . . ,  N}.  In  the  example  of 
Fig.  2,  n  =  5  and  the  children  are  Nodes  (2,  4,  5,  8,  9).  Using 
the  ray  from  Node  s  to  its  neighbor  with  lowest  ID  (Node  2) 
to  define  the  reference  direction,9  we  measure  the  angles 
{$2,4,  #2,5,  $2,8;  and  02,9}  (all  in  the  same,  e.g.,  counter-clock¬ 
wise,  direction),  as  shown  in  Fig.  2b.  The  complexity  of  this 
operation  is  0(n)  since  Node  s  has  n  children.  We  then  sort 
these  angles  in  increasing  order,  which  in  this  example  is: 
{$2,9,  #2, 4,  $2,8 ,  $2,5 }/  thereby  determining  that  the  "cyclic" 
ordering  of  the  nodes  is:  {2,  9,  4,  8,  5,  2).  This  sorting  process 
can  be  done  with  complexity  0(n  log  n)  using  heap  sort 
[16].  Once  this  is  done,  we  can  determine  the  nonoverlap¬ 
ping  angles  that  define  the  sectors  between  adjacent  nodes, 
namely,  (02.9,09,4,04,8,08,5,  and  05,2}-  We  then  find  the 
largest  such  angle,  which  in  this  case  is  84$  (resulting  in 
an  angle  of  coverage  of  0  =  360  —  04,g).  This  operation  has 
complexity  O(n).  Therefore,  the  overall  complexity  of  the 
sector-choosing  algorithm  for  a  node  with  n  children  (once 
the  node's  children  have  been  selected)  is 

0{n)  +  0(n  log  n)  +  O(n)  =  0(n  log  n). 

Finally,  since  the  number  of  children  a  node  may  have  is  at 
most  N  —  1,  the  overall  complexity  is  0(N  log  N)  at  each 
transmitting  node.  Since  there  are  O(N)  transmitting  nodes, 
the  beam  reduction  operation  (which  is  done  once  at  each 
transmitting  node  after  the  tree  has  been  constructed,  e.g., 
by  using  BIP)  has  overall  complexity  0(N 2  log  N).  Since 
BIP  has  complexity  0(N 3)  (see  Section  5.3),  adding  the 
beam-reduction  operation  (as  a  postalgorithm  process)  does 
not  affect  its  overall  degree  of  complexity,  and  the 
complexity  of  RB-BIP  remains  0(N3). 

5.2  An  Approach  Based  on  Incremental  Power: 
Directional  BIP  (D-BIP)  and  MIP  (D-MIP) 

In  [1]  and  [2],  we  proposed  the  Broadcast  Incremental 
Power  (BIP)  algorithm,  a  centralized  heuristic  for  the 
determination  of  low-power  broadcast  trees  in  networks 
with  omnidirectional  antennas.  BIP  is  the  basis  for  the 

9.  Any  child  node  could  have  been  used  to  determine  the  reference 
direction. 


Multicast  Incremental  Power  (MIP)  algorithm,  under  which 
the  tree  produced  by  BIP  is  pruned  by  eliminating  all 
transmissions  that  are  not  needed  to  reach  the  members  of 
the  multicast  group.  More  specifically,  under  MIP,  nodes 
with  no  downstream  destinations  do  not  transmit,  and 
some  nodes  may  be  able  to  reduce  their  transmitted  power 
(i.e.,  if  their  more-distant  downstream  neighbors  have  been 
pruned  from  the  tree). 

BIP  is  similar  in  principle  to  Prim's  algorithm  for  the 
formation  of  minimum-cost  spanning  trees  (MSTs),  in  the 
sense  that  new  nodes  are  added  to  the  tree  one  at  a  time  (on 
a  minimum-cost  basis)  until  all  nodes  are  included  in  the 
tree.  In  fact,  the  implementation  of  this  algorithm  is  based 
on  the  standard  Prim  algorithm,  with  one  fundamental 
difference.  Whereas  the  inputs  to  Prim's  algorithm  are  the 
link  costs  p,:l  (which  remain  unchanged  throughout  the 
execution  of  the  algorithm),  BIP  must  dynamically  update 
the  costs  at  each  step  (i.e.,  whenever  a  new  node  is  added  to 
the  tree).  This  updating  is  done  to  reflect  the  fact  that  the 
cost  of  adding  a  new  node  to  a  transmitting  node's  list  of 
neighbors  is  the  incremental  cost,  i.e.,  the  additional  cost 
associated  with  adding  a  new  downstream  neighbor,  given 
that  the  node  is  already  transmitting  at  some  particular 
power  level.  Consider  an  example  in  which  Node  i  is 
already  in  the  tree  (it  may  be  either  a  transmitting  node  or  a 
leaf  node)  and  Node  j  is  not  yet  in  the  tree.  If  Node  j  is 
already  participating  in  T  sessions  (hence,  no  transceivers 
are  available  for  an  additional  session),  the  cost  of  adding  it 
to  the  tree  is  set  to  00. 10  Otherwise,  for  all  such  Nodes  i  (i.e., 
all  nodes  already  in  the  tree)  and  Nodes  j  (i.e.,  nodes  not  yet 
in  the  tree),  the  following  is  evaluated: 

Pij  =  Pij~Pi,  (8) 

where  p,:l  is  the  link-based  cost  (power)  of  a  transmission11 
between  Node  i  and  Node  j  (i.e.,  it  is  r“  +  pT),  and  pt  is 
Node  i's  transmission  cost  prior  to  the  addition  of  Node  j; 
(which  includes  pT  if  node  i  is  already  transmitting;  if  Node 
i  is  currently  a  leaf  node.  Pi  =  0).  The  quantity  pb  represents 
the  incremental  cost  associated  with  adding  Node  j  to  the  set 
of  nodes  to  which  Node  i  already  transmits.  The  pair  {i,  j} 
that  results  in  the  minimum  value  of  pT  is  selected,  i.e.. 
Node  i  transmits  at  a  power  level  sufficient  to  reach  Node  j. 

10.  It  is  also  possible  to  associate  a  higher  cost  with  nodes  that  have  low 
"residual  capacity"  (i.e.,  few  available  transceivers);  however,  we  do  not  do 
so  in  this  paper. 

11.  We  neglect  pR  in  this  cost  measure  because  it  is  the  same  for  all 
possible  Node  js.  However,  pR  is  included  when  energy  consumption  is 
evaluated. 
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(C)  (d) 

Fig.  3.  Addition  of  nodes  in  D-BIP.  (a)  Source  node  and  three  destination  nodes,  (b)  Add  Node  2  to  tree,  (c)  Shift  beam  and  increase  power  to  add 
Node  3.  (d)  Use  Node  2  as  relay  to  reach  Node  1 . 


Thus,  one  node  is  added  to  the  tree  at  every  step  of  the 
algorithm. 

The  incremental  power  philosophy,  originally  developed 
for  use  with  omnidirectional  antennas,  can  be  applied  to 
tree  construction  in  networks  with  directional  antennas  as 
well.  At  each  step  of  the  tree-construction  process,  a  single 
node  is  added,  as  above.  However,  whereas  the  only 
variable  involved  in  computing  the  cost  (and  incremental 
cost)  in  the  omnidirectional  case  was  the  transmitter  power, 
the  directional-antenna  case  involves  the  choice  of  beam- 
width  9  as  well.  Based  on  the  propagation  model  of  (2),  the 
required  RF  power  increases  in  proportion  to  the  a  power 
of  the  distance  to  the  farthest  downstream  neighbor  and 
linearly  with  9. 

Consider  a  situation  in  which  Node  i  is  already 
transmitting  to  several  other  nodes.  The  incremental  cost 
of  adding  Node  j  to  Node  i's  set  of  downlink  neighbors 
depends  on  the  relative  location  of  Node  j  with  respect  to 
the  region  already  included  in  Node  i's  antenna's  cone  of 
coverage.  For  example,  if  Node  j  is  located  within  the  angle 
of  Node  i's  beam,  it  suffices  to  increase  Node  i's  commu¬ 
nication  range,  without  changing  the  width  or  direction  of 
the  beam.12  On  the  other  hand,  if  Node  j  is  not  located 
within  the  angle  of  Node  i's  beam,  then  the  beam  must  be 
adjusted;  this  is  usually  done  by  increasing  9,  although  it  is 
sometimes  possible  to  simply  shift  the  beam  if  all  of  a 
node's  downstream  neighbors  are  located  within  a  cone  not 
greater  than  0m;n.  Thus,  to  add  a  new  node,  it  is  sometimes 

12.  It  is  also  necessary  to  examine  whether  Node  j  could  be  added  to  the 
tree  at  lower  cost  by  using  a  different  node  (e.g.,  one  of  Node  i's 
downstream  neighbors)  as  its  upstream  neighbor. 


sufficient  to  simply  increase  transmission  range,  it  is 
sometimes  sufficient  to  simply  shift  the  beam,  sometimes 
the  beam  has  to  be  made  wider,  and  sometimes  a 
combination  of  increased  communication  range  and  beam 
characterization  must  be  done.  Note  that  there  is  no 
incremental  cost  associated  with  shifting  a  beam  (while 
maintaining  the  same  angle  of  coverage). 

Fig.  3a  shows  a  simple  example  in  which  the  source  node 
S  wants  to  add  nodes  1,  2,  and  3  to  the  tree.  Node  2  is  the 
closest  to  S,  so  it  is  added  first  (as  in  standard  omnidirec¬ 
tional  BIP);  in  Fig.  3b,  we  show  a  beam  of  width  0nlin  that  is 
centered  on  the  ray  between  S  and  Node  2.  We  must  then 
decide  which  node  to  add  next  (Node  1  or  Node  3),  and 
which  node  (that  is  already  in  the  tree)  should  be  its  parent. 
In  this  example,  the  distance  between  S  and  Node  3  is  only 
slightly  greater  than  that  between  S  and  Node  2;  addition¬ 
ally,  the  beam  from  S  to  Node  2  can  be  shifted  to  include 
both  Node  2  and  Node  3,  without  increasing  its  beamwidth. 
Therefore,  Node  3  is  added  next  by  increasing  the 
communication  range  of  S  and  by  shifting  the  beam  (see 
Fig.  3c).  Finally,  Node  1  must  be  added  to  the  tree.  One 
possibility  is  for  S  to  simply  increase  its  beamwidth, 
resulting  in  a  situation  similar  to  that  in  Fig.  2  (it  would 
not  have  to  increase  its  range  because  its  distance  to  Node  1 
is  less  than  its  distance  to  Node  3).  However,  doing  so 
would  result  in  a  huge  increase  in  beamwidth  (we  noted 
above  that  RF  power  is  directly  proportional  to  beam- 
width).  In  our  example,  it  is  more  cost-effective  for  Node  2 
to  become  the  parent  of  Node  1  (using  a  beamwidth  of  9mm, 
as  shown  in  Fig.  3d).  Note  that  the  tree  produced  by  this 
algorithm  depends  on  both  0m;n  and  the  propagation 
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constant  a.  By  contrast,  the  tree  produced  by  RB-BIP  is 
independent  of  0m;n,  although  it  does  depend  on  a. 

When  applied  to  the  broadcasting  problem,  the  resulting 
scheme  is  called  Directional  BIP  (D-BIP).  When  applied  to 
the  multicasting  problem,  the  first  step  is  to  form  a 
broadcast  tree  using  D-BIP.  To  implement  Directional  MIP 
(D-MIP),  the  broadcast  tree  produced  by  D-BIP  is  pruned, 
as  discussed  at  the  beginning  of  this  section.  Note  that  when 
^min  =  360,  D-BIP,  RB-BIP,  and  BIP  are  identical. 

In  [1],  [2],  we  described  the  "sweep"  operation,  in  which 
redundant  transmissions,  as  well  as  transmissions  that  can 
be  reduced  in  power  without  compromising  the  tree 
structure,  are  discovered.  Use  of  the  sweep  was  shown  to 
provide  modest  improvement  in  energy  efficiency  in 
networks  that  use  omnidirectional  antennas.  The  sweep 
operation  can  be  performed  either  after  the  tree  is 
constructed  (an  approach  suitable  for  all  tree-construction 
algorithms)  or  at  each  step  of  the  tree-construction  algo¬ 
rithm  (an  approach  feasible  only  for  those  algorithms,  such 
as  BIP,  that  add  a  single  node  to  the  tree  at  every  step). 
Typically,  more  improvement  is  observed  when  the  former 
approach  is  used  (in  omnidirectional  applications)  because 
the  impact  of  each  change  resulting  from  the  sweep  affects  a 
greater  number  of  nodes. 

However,  a  fundamental  characteristic  of  the  sweep 
makes  it  unattractive  in  applications  with  directional 
antennas.  In  omnidirectional  applications  (and,  hence, 
when  the  RB  versions  of  tree-construction  algorithms  are 
used  as  well),  the  sweep  often  results  in  an  increase  in  the 
number  of  one  or  more  nodes'  downstream  neighbors. 
Consequently,  there  is  often  an  increase  in  the  required 
angle  of  coverage  needed  by  the  antenna  in  RB  algorithms 
and,  thus,  in  an  increase  in  that  node's  RF  power,  and 
possibly  in  the  overall  tree  power.  With  D-BIP  and  D-MIP, 
this  problem  does  not  arise;  however,  there  are  fewer 
situations  in  which  the  sweep  will  discover  improvements 
(because  widening  the  beam  will  often  result  in  increased 
power).  Therefore,  the  performance  results  presented  in  this 
paper  do  not  make  use  of  the  sweep. 

5.3  Complexity  of  BIP  and  D-BIP 

First,  let  us  consider  the  complexity  of  BIP  with  omnidirec¬ 
tional  antennas.  One  node  is  added  to  the  tree  at  each  step; 
thus,  N  steps  are  needed  to  construct  a  tree  with  N  nodes. 
Consider  step  k,  at  which  k  nodes  are  already  in  the  tree  (we 
call  them  "inside"  nodes)  and  the  remaining  (N  —  k) 
"outside"  nodes  must  eventually  be  added.  Thus,  the 
number  of  inside  nodes  is  O(N),  and  the  number  of  outside 
nodes  is  also  O(N).  At  each  step,  the  algorithm  must 
determine  which  node  can  be  added  at  least  incremental 
cost.  Thus,  it  is  necessary  to  know  the  cost  from  each  inside 
node  to  every  outside  node,  i.e.,  for  each  inside  node,  we 
must  find  the  least-cost  outside  node;  we  then  take  the 
minimum  over  all  inside  nodes. 

Most  of  the  database  needed  to  implement  BIP  remains 
unchanged  from  one  step  to  the  next.  The  only  costs  that 
must  be  updated  are  those  that  relate  to  the  new  child  node 
(which  changes  its  status  from  an  outside  node  to  an  inside 
node)  and  to  its  parent.  The  cost  between  the  new  child 
node  and  all  of  the  remaining  outside  nodes  must  be 
computed;  the  complexity  of  this  operation  is  O(N). 


Additionally,  since  its  parent  has  increased  its  RF  power, 
the  cost  from  the  parent  to  all  remaining  outside  nodes 
must  be  updated;  this  also  has  complexity  O(N).  For  the 
other  inside  nodes,  we  simply  delete  the  entry  correspond¬ 
ing  to  the  node  that  was  just  added  to  the  tree;  this  also  has 
complexity  O(N).  Therefore,  updating  the  database  re¬ 
quires  three  sequential  operations  of  complexity  O(N), 
which  results  in  a  complexity  of  O(N).  Once  the  database  is 
updated,  the  minimum-cost  new  node  must  be  chosen.  The 
complexity  of  this  search  is  at  most  0(N2).  Therefore,  the 
complexity  at  each  step  is  at  most  0(N2).  Since  there  are  N 
steps,  the  overall  complexity  is  at  most  0(N3). 

Now,  let  us  consider  the  complexity  of  D-BIP,  for  which 
we  must  incorporate  the  impact  of  the  choice  of  antenna 
beam  angle  at  each  step  (in  contrast  to  the  case  of  RB-BIP,  in 
which  the  beam  reduction  process  was  executed  once  at 
each  transmitting  node,  after  the  tree  was  constructed).  As 
in  the  case  of  omnidirectional  BIP,  most  of  the  database 
remains  unchanged  from  one  step  to  the  next.  The 
recomputation  of  costs  to  outside  nodes  is  necessary  only 
for  two  nodes,  namely,  the  parent  of  the  node  that  was  just 
added  and  for  the  child  (since  it  becomes  an  inside  node). 
First,  consider  the  parent.  Following  the  discussion  of 
Section  5.1,  the  complexity  of  choosing  the  best  outside 
node  for  any  given  inside  node  (which  involves  choosing 
the  antenna  beamwidth  6  as  well)  is  0(N2  log  N).  Now, 
consider  the  child.  It  is  necessary  to  compute  its  cost  to  all 
outside  nodes.  Since  the  child  doesn't  have  any  children  of 
its  own  yet,  it  is  not  necessary  to  compute  the  beamwidth  6 
because  it  is  0m;n  whenever  a  node  has  a  single  child.  Thus, 
the  complexity  of  computing  the  cost  from  the  child  to  all 
remaining  outside  nodes  O(N).  For  every  other  inside  node, 
the  newly  added  node  is  simply  deleted  from  its  list  of 
outside  nodes,  without  recomputing  anything.  Therefore, 
the  complexity  of  each  step  is  0(N2  log  N)  and  since  there 
are  N  steps,  the  overall  complexity  of  D-BIP  is  0(N3  log  N). 

5.4  Example  Broadcast  Trees 

Fig.  4a  shows  the  broadcast  tree  produced  by  BIP  for  a 
10-node  network,  where  the  source  node  is  shown  larger 
than  the  other  nodes.  As  noted  in  Section  3.1,  RB-BIP  uses 
the  same  tree  as  BIP  (which  is  based  on  omnidirectional 
antennas);  the  only  difference  is  that  the  antenna  beam- 
width  is  reduced.  Fig.  4b  shows  the  optimal  tree  for 
omnidirectional  antennas,  which  was  obtained  by  exhaus¬ 
tive  search.  The  tree  structure,  as  well  as  the  resulting  value 
of  total  tree  power  P,  depend  on  the  value  of  the 
propagation  constant  a;  our  results  are  based  on  a  =  2. 
Tree  power  P  is  listed  in  the  figure  caption  for  0niin  =  360 
(the  omnidirectional  case),  as  well  as  0niin  =  30  and  1.  There 
is  relatively  little  power  savings  when  6>nljn  is  reduced  below 
30  because  the  two  highest-range  transmissions  require  the 
use  of  6  >  30  to  reach  all  of  their  downstream  neighbors. 

Under  D-BIP  (unlike  RB-BIP),  the  tree  structure  depends 
on  the  value  of  0,mn.  Fig.  5a  shows  the  tree  for  the  same 
network  for  D-BIP  with  0„un  =  1.  In  this  example,  D-BIP 
produces  a  tree  in  which  each  node  has  only  a  single 
downstream  neighbor  (thus,  6  =  9min  at  each  node)  result¬ 
ing  in  a  zigzag  path  with  no  branching.  The  value  of  P  is 
greatly  reduced  by  using  highly  directional  antennas. 
However,  this  value  is  84  percent  greater  than  that  of  the 
optimal  tree  for  f?m;n  =  1,  as  shown  in  Fig.  5b. 
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Fig.  4.  Example  1 0-node  broadcast  trees  based  on  use  of  omnidirectional  antennas  (the  same  tree  is  used  for  RB-BIP,  independent  of  the  value  of  0niin). 
(a)  BIP,  RB-BIP:  0mill  =  360  :  P  =  14.06,  6Un  =  30  :  P  =  4.26,  and  6»mill  =  1  :  P  =  3.99.  (b)  optimal  (based  on  6»min  =  360).  6»min  =  360  :  P  =  10.71, 
6U„  =  30  :  P  =  3.728,  and  6U,  =  1  :  P  =  3.709. 


Fig.  5.  Example  10-node  broadcast  trees  for  9^  =  1.  (a)  D-BIP  (P  =  0.1051)  and  (b)  optimal  ( P  =  0.05707). 


Fig.  6.  Example  10-node  broadcast  trees  for  9^  =  30.  (a)  D-BIP  ( P  =  1.722)  and  (b)  optimal  ( P  =  1.607). 


Fig.  6a  shows  the  D-BIP  tree  for  the  same  network,  but 
with  0m in  =  30.  Here,  P  is  7.2  percent  greater  than  that  of  the 
optimal  tree,  as  shown  in  Fig.  6b. 

These  results  demonstrate  that  the  use  of  directional 
antennas  can  facilitate  considerable  energy  saving  through 
the  use  of  algorithms  such  as  RB-BIP  and  D-BIP.  More¬ 
over,  D-BIP  provides  lower-power  trees  than  RB-BIP  for  a 
given  value  of  9m ;n,  and  this  advantage  increases  as  6m\n 
decreases.  However,  when  7>rnill  is  very  small,  even  the  tree 
produced  by  D-BIP  is  likely  to  have  a  significantly  higher 
value  of  RF  transmission  power  than  the  optimal  tree  (on 
a  percentage  basis). 


We  attribute  the  relatively  good  performance  of  BIP 
when  0nljn  >  6  (as  measured  by  the  closeness  of  tree  power 
to  its  optimal  value,  on  a  percentage  basis)  to  the  wireless 
multicast  advantage  (see  Section  3.1).  However,  this 
property  no  longer  applies  when  highly  directional  anten¬ 
nas  are  used  because  power  is  directly  proportional  to 
beamwidth  9 ;  thus,  it  is  costly  to  expand  a  beam  to 
accommodate  additional  nodes.  Therefore,  the  greedy 
nature  of  our  incremental  power  approach  suffers  when 
used  with  extremely  narrow  beams,  and  alternative 
approaches  may  be  desirable. 
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Fig.  7.  Average  RF  tree  power  for  broadcasting  examples  with  a  50- 
node  network  using  BIP,  based  on  use  of  directional  antennas 

( Pmax  =  25,  a  =  2). 

5.5  Average  Tree  Power 

We  have  simulated  BIP  and  MIP  using  the  two  directional 
antenna  schemes;  thus,  we  have  simulated  a  total  of  four 
algorithms,  namely,  RB-BIP,  D-BIP,  RB-MIP,  and  D-MIP.  In 
all  cases,  we  consider  a  network  of  N  =  50  nodes  that  are 
randomly  located  in  a  region  with  dimensions  5x5 
(arbitrary  units  of  distance);  the  same  node  locations  are 
used  in  our  broadcasting  and  multicasting  examples.  We  let 
a  =  2,  i.e.,  the  required  RF  power  value  is  r2.  RF  transmis¬ 
sion  power  levels  are  bounded  by  pmin  =  0  and  pmax  =  25 
(corresponding  to  a  maximum  communication  range  of  5). 
The  results  presented  in  this  section  are  based  on  the 
availability  of  unlimited  numbers  of  transceivers,  frequen¬ 
cies,  and  energy;  thus,  it  is  always  possible  to  include  all 
destinations  in  the  tree.  We  are  thus  able  to  focus  on  the 
impact  of  adaptive  beamwidth  capabilities  on  energy 
efficiency.  We  also  neglect  the  impact  of  signal-processing 
power  by  setting  (pT,pR)  =  (0,0).  Our  results  are  based  on 
simulation  runs,  each  consisting  of  X  =  10, 000  multicast 
sessions. 

First,  we  consider  broadcasting.  The  source  node  is 
chosen  at  random  for  each  session,  and  RB-BIP  or  D-BIP 
is  used  to  construct  the  broadcast  tree  based  on  that 
source  node.  Fig.  7  shows  the  average  tree  RF  power  as  a 
function  of  0mm  for  RB-BIP  and  D-BIP.  As  expected,  both 
schemes  provide  substantial  reduction  in  RF  power  as 
compared  to  omnidirectional  antennas  when  0mm  360 
because  energy  is  concentrated  where  it  is  needed.  For 
Omin  >  180  there  is  vitually  no  difference  in  performance 
between  the  two  schemas.  However,  for  small  values  of 
f?min,  D-BIP  provides  a  substantial  reduction  in  RF  power 
as  compared  with  RB-BIP. 

Fig.  8  shows  average  tree  power  for  multicasting 
examples  in  which  RB-MIP  and  D-MIP  are  used.  As  above, 
the  source  node  is  chosen  at  random  for  each  session. 
Multicast  groups  are  chosen  randomly  for  each  session 
request;  the  number  of  destinations  is  uniformly  distributed 
between  1  and  N  —  1.  Results  are  qualitatively  similar  to 
those  for  the  broadcasting  case,  although  the  values  of  tree 
power  are  somewhat  lower.  The  power  of  a  multicast  tree 
generated  by  MIP  is  always  less  than  or  equal  to  that  of  a 


Fig.  8.  Average  RF  tree  power  for  multicasting  examples  with  a  50-node 
network  using  MIP,  based  on  use  of  directional  antennas 

(Pmax  =  25,  a  =  2). 

broadcast  tree  generated  by  BIP  because  MIP  forms  multi¬ 
cast  trees  by  pruning  broadcast  trees  generated  by  BIP. 

In  view  of  the  similarity  (both  qualitative  and  quantita¬ 
tive)  of  the  results  for  broadcasting  and  multicasting,  we 
focus  on  multicasting  in  the  rest  of  this  paper. 

6  The  Incorporation  of  Resource  Limitations 

The  discussions  in  the  previous  sections  implicitly  assume 
that  sufficient  resources  are  available  to  implement  the  trees 
created  by  the  algorithms.  These  resources  include  trans¬ 
ceivers,  frequencies,  and  battery  energy.  In  this  section,  we 
discuss  how  limitations  on  these  resources  are  incorporated 
into  our  model  and  how  our  algorithms  can  be  modified  to 
cope  with  limited  energy. 

It  is  straightforward  to  incorporate  the  impact  of  a  finite 
number  of  transceivers.  When  constructing  a  tree  for  a  new 
arrival,  the  cost  of  a  node  is  set  to  oc  if  all  of  its  transceivers 
are  currently  supporting  other  sessions.  However,  the 
modeling  of  finite  frequency  resources  is  much  more 
complicated. 

6.1  Bandwidth  Limitations 

Let  us  consider  the  case  in  which  Node  m  wants  to  transmit 
to  Node  n.  Any  particular  frequency /may  be  unusable  for 
one  of  the  following  reasons: 

•  /is  already  in  use  (for  either  transmission  or 
reception)  at  either  Node  m  or  Node  n, 

•  /is  being  used  by  one  or  more  nodes  that  create 
interference  at  Node  n,  thereby  preventing  reception 
at/,  and 

•  the  use  of/by  Node  m  would  interfere  with  ongoing 
communications  at  other  nodes. 

In  this  paper,  we  use  the  following  basic  greedy 
approach  for  frequency  assignment,  which  we  referred  to 
as  FA1  in  [13]: 

Assume  the  availability  of  an  infinite  number  of 
frequencies  when  forming  the  tree  (the  approach  used  in 
[1]  and  [14]).  Then  attempt  to  assign  the  available 
frequencies  to  the  tree.  The  assignment  process  is  complete 
when  either  frequencies  have  been  assigned  to  all  transmis- 
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sions  or  when  no  additional  frequencies  are  available  to 
support  portions  of  the  tree. 

Under  this  scheme,  the  tree  construction  process  ignores 
the  possibility  that  frequencies  may  not  be  available  to 
provide  the  required  connectivity.  Thus,  if  appropriate 
frequencies  cannot  be  found  along  the  paths  to  all  desired 
destinations,  then  some  destinations  will  not  be  reached. 
We  have  used  a  greedy  version,  in  which  frequencies  are 
assigned  using  an  orderly  procedure,  without  the  possibi¬ 
lity  of  backtracking  to  change  assignments  and  without  the 
use  of  exhaustive  search  (or  other  scheme)  to  determine 
whether  a  consistent  frequency  assignment  is  possible. 
Specifically,  we  simply  assign  the  lowest-numbered  avail¬ 
able  noninterfering  frequency  to  each  node.  Thus,  this 
scheme  can  result  in  unreached  destinations,  even  though 
they  might  be  reachable  through  a  better  frequency  assign¬ 
ment.  But,  this  is  a  common  characteristic  of  all  heuristic 
procedures. 

In  [13],  we  also  considered  an  alternative  scheme  (FA2) 
under  which,  at  each  step  of  the  tree-construction,  the 
frequency  is  chosen  along  with  the  transmission  power 
level.  Under  FA2,  the  tree  is  formed  using  only  nodes  that 
do,  in  fact,  have  frequencies  available.  Again,  there  is  no 
guarantee  that  all  destinations  will  be  reached.  However, 
FA2  provides  a  richer  search  space  than  FA1.  In  this  paper, 
we  focus  exclusively  on  FA1  because  it  is  simple  to  use  and 
is  applicable  to  any  tree-construction  algorithm.  FA2  can  be 
used  with  BIP  (and  similar  schemes  in  which  one  node  is 
added  to  the  tree  at  each  step),  but  not  with  some  of  the 
other  algorithms  discussed  in  [1]  and  [2]. 


6.2  Energy  Limitations 

Use  of  a  cost  metric  that  involves  only  the  total  power 
required  to  maintain  the  tree  can  result  in  rapid  energy 
depletion  at  some  nodes.  When  nodes  "die"  in  this  manner, 
it  may  be  no  longer  possible  to  create  energy-efficient  trees. 

We  can  discourage  the  inclusion  of  energy-poor  nodes  in 
the  multicast  tree  by  increasing  the  cost  associated  with 
their  use.  In  (4),  we  defined  the  residual  energy  at  Node  i  at 
time  t  to  be  We  now  define  the  cost  of  a  link  between 
Node  i  and  Node  j  to  be 


Cij  —  Pij 


EMV 


(9) 


where  (3  is  a  parameter  that  reflects  the  importance  we 
assign  to  the  impact  of  residual  energy.13  Clearly,  when 
(3  =  0,  the  link  cost  is  simply  the  power  needed  to  maintain 
the  link. 

The  incremental  cost  associated  with  adding  Node  j  to 
the  set  of  Node  i's  downstream  neighbors,  given  that  Node  i 
is  already  transmitting  at  power  level  pt  (hence,  at  cost  C() 
is: 


when  (3  is  too  large,  too  much  emphasis  may  be  placed  on 
balancing  energy  use  throughout  the  network,  while 
underemphasizing  the  need  for  energy  efficiency. 

Performance  results  in  Section  7  and  [8]  show  the 
beneficial  effects  of  using  (3  in  the  range  [0.5,  2].  It  would  be 
possible  to  develop  alternative  cost  functions  to  (9)  that  also 
discourage  the  use  of  energy-poor  nodes;  we  make  no  claim 
of  optimality.  Our  objective  is  to  demonstrate  that  load 
balancing  based  on  residual  energy  can  extend  a  network's 
useful  lifetime. 

7  Performance  Results 

Important  performance  measures  for  energy-constrained 
networks  include  network  lifetime  and  delivered  traffic 
volume.  In  this  section,  we  present  our  performance  results 
for  the  two  multicast  tree-construction  algorithms  we  have 
developed  for  directional  antennas,  namely,  Reduced- 
Beamwidth  MIP  (RB-MIP)  and  Directional  MIP  (D-MIP). 

As  in  Section  5.5,  we  have  simulated  the  performance  of 
RB-MIP  and  D-MIP  for  a  network  of  N  =  50  nodes  that  are 
randomly  located  in  a  region  with  dimensions  5x5;  the 
same  node  locations  are  used  in  all  examples  presented  in 
this  paper.  In  extensive  performance  evaluation,  we  have 
observed  that  these  results  are  representative  of  other 
random  node  distributions  as  well.  We  present  results 
primarily  for  a  propagation  constant  value  of  a  =  2,  which 
results  in  required  RF  power  values  of  r2  to  support  a  link 
between  two  nodes  that  are  separated  by  distance  r;  we  also 
present  a  limited  set  of  results  for  a  =  4.  We  set  arbitrary 
values  for  transmission  processing  power  (pT )  and  recep¬ 
tion  processing  power  (pR ).  In  particular,  we  consider 
(pT,pR)  =  (0,0)  as  well  as  "moderate"  (0.01,  0.1)  and  "high" 
(0.1,  1)  values  of  these  quantities.  RF  transmission  power 
levels  are  bounded  by  pmin  =  0  and  pmax  =  25.  In  most  of  our 
experiments,  the  initial  energy  at  each  node  is  200  (arbitrary 
units,  consistent  with  the  units  of  distance).14  We  demon¬ 
strate  the  impact  of  incorporating  residual  energy  into  the 
cost  metric  and  compare  performance  for  (3=0,  0.5,  1,  and 
2. 

In  our  simulations,  multicast  requests  arrive  with 
interarrival  times  that  are  exponentially  distributed  with 
rate  X/N  at  each  node;  we  have  used  A  =  1  in  our 
simulations.  Session  durations  are  exponentially  distributed 
with  mean  1.  Multicast  groups  are  chosen  randomly  for 
each  session  request;  the  number  of  destinations  is  uni¬ 
formly  distributed  between  1  and  N  —  1. 

Each  simulation  run  consists  of  X  =  10,  000  multicast 
sessions,  some  of  which  may  be  blocked  because  of  lack  of 
resources  (which  in  general  include  transceivers,  frequen¬ 
cies,  and  energy).  The  same  random  number  sequence  is 
used  to  drive  each  of  our  experiments,  thereby  facilitating  a 
meaningful  comparison  of  results  for  different  values  of  (3. 


c'j-Qj  a.  (io) 

When  (3  is  too  small,  too  much  emphasis  may  be  placed 
on  the  construction  of  energy  efficient  trees,  resulting  in  the 
rapid  depletion  of  energy  at  some  of  the  nodes.  By  contrast, 

13.  Residual  energy  was  incorporated  into  the  cost  metric  in  a  similar 
manner  in  [17]. 


7.1  Network  Lifetime 

A  fundamental  issue  in  limited-energy  applications  is  net¬ 
work  lifetime,  i.e.,  the  interval  over  which  the  network  can 

14.  We  assume  that  if  a  node  is  alive  at  the  beginning  of  a  session,  it  will 
be  able  to  complete  the  session  (regardless  of  whether  it  is  a  transmitting  or 
a  receive /only  node).  Thus,  we  neglect  the  minor  "end  effects"  associated 
with  a  node's  death  during  a  session. 
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Fig.  9.  Evolution  of  number  of  live  nodes  under  MIP  with  omnidirectional 
antennas  for  50-node  network  (zero  processing  power). 

provide  acceptable  levels  of  service.  Clearly,  a  suitable 
definition  of  network  lifetime  depends  on  the  specific 
application.  For  example,  in  some  applications,  one  may 
view  network  death  as  the  time  at  which  the  first  node  dies 
(e.g.,  see  [17])  because  it  is  no  longer  possible  to  reach  all  of  the 
nodes.  Alternatively,  network  death  may  be  defined  as  the 
death  of  a  specified  fraction  of  the  nodes.  In  this  paper,  we 
don't  specify  a  particular  definition  of  network  death, 
although  we  do  feel  that  a  reasonable  definition  of  acceptable 
performance  would  require  that  at  least  50  percent  of  the 
nodes  remain  alive.  Instead,  we  examine  the  time  evolution  of 
the  number  of  live  nodes. 

In  this  section,  we  consider  the  case  of  unlimited 
numbers  of  transceivers  and  frequencies,  but  finite  energy 
at  each  node.  We  present  results  for  omnidirectional 
antennas  (although  results  are  qualitatively  similar  for 
directional  antennas).  Thus,  we  are  able  to  focus  on  the 
impact  of  energy  constraints,  without  addressing  other 
system  parameters.  In  such  cases,  all  desired  destinations 
can  be  reached,  provided  that  live  nodes  are  available  to 
support  the  required  trees. 

Fig.  9  shows  the  evolution  of  the  number  of  live  nodes  as  a 
function  of  the  number  of  session  arrivals  for  (3  =  0,  0.5,  1, 
and  2.  Results  are  shown  for  the  case  of  zero  processing 
power,  i.e.,  (pT ,pR)  =  (0, 0).  As  noted  in  Section  6,  the  use  of 
nonzero  values  of  / 3  tends  to  discourage  the  use  of  nodes  that 
have  little  residual  energy.  The  use  of  0.5  <  (3  <  2,  rather 
than  0,  results  in  a  significant  delaying  of  the  first  node's 
death,  and  keeps  a  large  fraction  (e.g.  80  percent  or 
90  percent)  of  the  nodes  alive  for  a  considerably  greater 
number  of  sessions.  Specifically,  for  zero  precessing  power, 
when  (3  =  0,  the  first  node  dies  at  arrival  136;  for  f3  =  0.5, 1, 
and  2,  the  first  node  dies  at  arrival  563,  668,  and  599, 
respectively.  Results  are  qualitatively  similar  when 
( pT,pR )  =  (0.1, 1),  except  that  nodes  die  much  faster  because 
of  the  energy  consumed  by  signal  processing  (see  [8]). 

Moreover,  for  0.5  <  (3  <  2,  once  about  10  percent  of  the 
nodes  have  died,  the  fraction  of  live  nodes  decreases  to 
below  10  percent  shortly  thereafter.  The  rapid  death  of 
nodes  in  this  manner  is  not  a  harmful  effect.  Once  about 
50  percent  of  the  nodes  are  dead,  a  significant  number  of  the 
remaining  live  nodes  are  typically  unreachable.  Thus,  the 
fact  that  the  use  of  (3  =  0  maintains  a  certain  fraction  (say 
25  percent)  of  the  nodes  alive  considerably  longer  than  use 
of  larger  values  of  (3  is  not  beneficial. 


Fig.  10.  Evolution  of  cumulative  bit  volume  under  MIP  for  several  sets  of 
(F,  T)  pairs  (/?=  0 ;(pr,pR)  =  (0,0)). 

Thus,  for  0.5  <  (3  <  2,  we  have  achieved  a  high  degree  of 
load  balancing  that  keeps  almost  all  of  the  nodes  alive  for  a 
relatively  long  time,  thereby  maintaining  network  connec¬ 
tivity  and  high  levels  of  throughput  much  longer  than  for 
the  case  in  which  f3  =  0.  In  view  of  the  relative  insensitivity 
of  node  lifetime  to  the  value  of  (3  (in  the  region  0.5  <f3<  2), 
we  use  (3  =  1  in  the  examples  presented  in  this  paper.  No 
claim  for  optimality  is  made. 

Since  we  use  a  finite  value  of  pmax,  it  is  typical  to  achieve 
a  final  state  in  which  a  number  of  nodes  still  have  energy, 
but  further  communication  is  impossible  because  of  a  lack 
of  connectivity  among  the  live  nodes.15 

7.2  Delivered  Traffic  Volume 

We  now  consider  the  delivered  traffic  volume  Btotal.  In 
doing  so,  we  address  the  impact  of  realistic  constraints  on 
the  number  of  transceivers  (T)  available  at  each  node  and  on 
the  number  of  frequencies  (F)  available  for  communication. 
Our  modeling  assumptions  are  the  same  as  those  of  the 
previous  section.  Unlike  the  case  of  infinite  transceiver  and 
frequency  resources,  performance  depends  strongly  on  the 
arrival  rate  A  because  high  traffic  loads  require  a  large 
number  of  transceivers  and  frequencies  to  support  them. 
We  present  results  for  MIP,  first  for  omnidirectional  and 
then  for  directional  antennas.  Our  results  are  based  on  the 
use  of  frequency  assignment  scheme  FA1.  As  noted  in 
Section  5.2,  the  "sweep"  is  not  used. 

7.2. 1  Omnidirectional  Antennas 

Fig.  10  shows  the  time  evolution  of  Btotal  under  MIP,  with 
omnidirectional  antennas.  To  illustrate  the  impact  of  a  finite 
number  of  transceivers  and  frequencies,  we  present  results 
for  ( F,T )  =  (4,2)  and  (oo,oo).  To  illustrate  the  impact  of 
using  a  nonzero  value  of  (3,  we  present  results  for  (3  =  0  and 
1.  These  results  are  based  on  an  arrival  rate  of  A  =  1  and 
(pT ,pR)  =  (0, 0).  One  unit  on  the  vertical  axis  corresponds  to 
the  delivery  of  a  message  of  average  length  (one  time  unit) 
to  a  single  destination  (see  definition  in  (5)).  The  initial 
value  of  energy  at  each  node  is  £)(0)  =  200. 

Consider  first  the  curves  for  (3  =  0.  The  use  of  (F,  T)  = 
(4, 2)  results  in  a  significantly  slower  rise  in  the  value  of 
delivered  traffic  volume;  the  blocking  of  destinations  (or 

15.  If  we  had  included  the  energy  consumption  associated  with 
operation  when  a  node  is  simply  "on,"  but  neither  transmitting  nor 
receiving,  all  nodes  would  have  died  eventually. 
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number  of  arrivals 

Fig.  11.  Evolution  of  cumulative  bit  volume  under  MIP  with 
directional  antennas  for  D-MIP  and  RB-MIP;  T  =  oo,  F  =  oo 

(/3=l;(pT,pfl)  =  (0,0)). 

perhaps  even  entire  sessions)  results  from  the  unavailability 
of  transceiver  or  frequency  resources  at  one  or  more  nodes. 
Nevertheless,  the  final  value  (i.e.,  the  total  delivered  traffic 
volume)  is  4.16  percent  greater  than  that  for  the  case  of 
( F,T )  =  (oo,  oo).  Overall,  the  value  of  Btotal  is  relatively 
insensitive  to  the  values  of  F  and  T,  even  though  these 
parameters  have  a  significant  impact  on  the  rate  of  traffic 
delivery.  A  possible  explanation  for  the  higher  value  of 
Btotal  js  achieved  for  low  values  of  F  and  T  is  that  costly 
destinations  are  more  likely  to  be  blocked  when  few 
resources  are  available,  thus  resulting  in  a  lowered  average 
cost  per  destination  that  is  actually  reached. 

Now,  let  us  consider  the  impact  of  setting  (3  =  1.  For  any 
(F,T)  pair,  the  curve  can  be  approximated  well  by  a  linear 
increase  until  the  final  value  is  reached,  a  departure  from 
the  asymptotic  performance  observed  for  (3  =  0.  This 
behavior  can  be  explained  by  the  fact  that  the  load 
balancing  that  results  from  the  use  of  (3  =  1  provides  a 
rapid  transition  from  a  state  in  which  most  nodes  are  alive 
to  one  in  which  most  are  dead,  as  shown  in  Fig.  9.  Thus, 
there  are  two  distinct  regions  of  operation.  When  all  (or 
most)  nodes  are  alive,  the  rate  of  traffic  delivery  is 
maintained  at  (or  near)  its  maximum  value.  When  most 
nodes  are  dead,  the  rate  of  traffic  delivery  is  close  to  (or 
equal  to)  zero.  We  again  observe  that,  although  ( F,T )  = 
(4, 2)  provides  a  lower  initial  rate  of  traffic  delivery,  it 
provides  an  8.22  percent  higher  value  of  Btotal.  Additional 
discussion  of  delivered  traffic  volume  in  energy-con- 
strained  systems  with  omnidirectional  antennas  may  be 
found  in  [8]. 

7.2.2  Directional  Antennas 

We  now  consider  the  case  of  directional  antennas.  Fig.  11 
shows  the  time  evolution  of  Btotal  for  RB-MIP  and  D-MIP  for 
several  values  of  0mm.  Results  are  shown  for  (3=1 ,  zero 
processing  power,  and  T  =  F  =  oo.  The  case  of  0m;n  =  360 
corresponds  to  the  use  of  omnidirectional  antennas.  Our 
first  observation  is  that  the  use  of  RB-MIP  and  D-MIP 
provide  significantly  increased  values  of  delivered  traffic 
volume,  and  that  this  volume  increases  as  0m;n  decreases. 
The  increase  is  less  than  linear  in  l/#min  because  some 
beamwidths  may  be  greater  than  6 

For  0m;n  =  30,  60,  and  90,  two  curves  are  shown  for  each 
value;  the  lower  curve  is  for  RB-MIP  and  the  upper  curve  is 


Fig.  12.  Evolution  of  cumulative  bit  volume  under  MIP  with  directional 
antennas  for  D-MIP  and  RB-MIP;  T  =  2,  F  =  4  (/3  =  1;  (jpr,pR)  =  (0,0)). 

for  D-MIP.  In  all  cases,  D-MIP  provides  better  performance 
than  RB-MIP,  and  its  advantage  increases  as  0n,m  decreases. 
As  in  the  omnidirectional  case  for  (3=1 ,  the  curves  can  be 
closely  approximated  by  straight  lines  with  slope  indepen¬ 
dent  of  $mm,  until  the  final  value  is  reached.  For  the  present 
case  of  unlimited  transceiver  and  frequency  resources,  this 
observation  is  not  surprising.  Since  the  network  has 
sufficient  resources  to  deliver  all  offered  traffic,  and  since 
the  use  of  (3  =  1  achieves  good  load  balancing  (causing  the 
nodes  to  die  at  almost  the  same  time),  the  delivered  traffic 
volume  increases  essentially  linearly  until  the  nodes  die. 
Reducing  the  value  of  f?m;n  results  in  a  reducing  of  the 
energy  expenditure,  while  not  affecting  the  rate  of  traffic 
delivery  (all  offered  traffic  while  the  nodes  are  alive,  and 
none  while  they  are  dead).  However,  the  point  at  which  the 
nodes  die  (as  a  function  of  9mm)  cannot  be  predicted  without 
actually  running  the  simulation. 

Fig.  12  shows  similar  results  for  finite  transceiver  and 
frequency  resources,  namely  ( F,T )  =  (4,2).  As  in  the  case 
of  infinite  resources,  each  curve  shows  a  linear  increase 
until  its  final  value  is  reached;  however,  there  are  noticeable 
differences  in  slope  for  different  values  of  9mi„-  As  a 
consequence  of  the  reduced  resources,  a  significant  fraction 
of  destinations  is  blocked,  resulting  in  a  slower  approach 
(lower  slope)  to  the  final  value  of  total  delivered  traffic 
volume.  Nevertheless,  the  final  value  differs  little  from  that 
for  the  case  of  infinite  transceiver  and  frequency  resources, 
as  we  observed  for  the  case  of  omnidirectional  antennas. 

Table  1  shows  the  total  delivered  traffic  volume 
(corresponding  to  the  final  values  of  the  curves  in  Figs.  11 
and  12)  for  both  RB-MIP  and  D-MIP,  four  values  of  0mm, 
and  three  (F,T)  pairs.  This  table  shows  that  the  impact  of 
finite  resources  on  total  delivered  traffic  volume  is  small, 
even  though  the  rate  at  which  traffic  is  delivered  (the  slope 
of  the  curve)  is  reduced  significantly. 

Fig.  13  shows  the  impact  of  a  higher  value  of  the 
propagation  constant,  namely  a  =  4,  for  unlimited  trans¬ 
ceiver  and  frequency  resources.  It  is  appropriate  to  consider 
higher  values  of  a  such  as  this  because  the  propagation  loss 
rate  increases  as  the  transmitting  antenna  is  moved  closer  to 
ground  level  and  as  the  density  of  foliage  increases.  The 
curves  are  qualitatively  similar  to  those  of  Fig.  11,  in  that  the 
curves  can  again  be  closely  approximated  by  straight  lines 
until  the  final  value  is  reached,  again  with  slope  that  is 
independent  of  9m;n.  As  in  the  case  for  a  =2,  there  is  a 
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TABLE  1 

Total  Delivered  Traffic  Volume;  a  =  2 


emin  =  360 

0mill  =  90 

0mi„  =  60 

emill  =  30 

(F.  T ) 

RB-MIP 

D-MIP 

RB-MIP 

D-MIP 

RB-MIP 

D-MIP 

(oo;  oo) 

19.447 

64,396 

68.732 

79,287 

95,541 

103,518 

146,821 

(8.4) 

20,364 

64,805 

68,945 

79.118 

95.789 

102,997 

146,761 

(4.  2) 

21.046 

64,380 

71,169 

78,522 

96,129 

100,625 

150,347 

strong  dependence  of  the  final  value  of  delivered  traffic 
volume  on  However,  there  is  significantly  less 

difference  between  the  results  for  D-MIP  and  RB-MIP.  For 
example,  for  0nlin  =  30,  this  difference  is  only  16.97  percent, 
as  compared  to  41.83  percent  for  the  corresponding  case  of 
a  =  2  and  unlimited  resources. 

The  impact  of  increasing  a  can  be  explained  qualitatively 
by  recognizing  that,  as  a  increases,  the  penalty  for  using 
longer  links  increases.  Thus,  trees  typically  consist  of  a 
larger  number  of  shorter  links.  Consequently,  under 
omnidirectinal  MIP  (and,  hence,  RB-MIP  as  well)  nodes 
will  tend  to  have  a  smaller  number  of  downstream 
neighbors  and,  hence,  the  average  beamwidth  used  for 
RB-MIP  will  tend  to  be  less  than  that  for  smaller  values  of  a. 
Therefore,  there  will  be  less  of  a  difference  in  the  trees 
produced  by  RB-MIP  and  D-MIP  (which  tends  to  produce 
trees  in  which  a  node's  downstream  neighbors  are  located 
within  relatively  narrow  beams). 

When  comparing  results  for  a  =  2  and  a  =  4,  we  do  not 
attach  any  significance  to  the  difference  in  total  delivered 
traffic  volume  or  to  network  lifetime.  It  is  not  possible  to 
compare  directly  the  results  for  different  values  of  a 
because  we  use  arbitrary  units  of  distance,  and  because 
different  constant  factors  may  be  associated  with  propaga¬ 
tion  loss  for  different  values  of  a  in  realistic  environments. 

7.3  Traffic  Volume  per  Unit  Energy 

It  is  also  of  interest  to  study  the  dependence  of  traffic 
volume  on  9min.  Fig.  14  shows  Bx,e  the  total  number  of  bits 
delivered  per  unit  energy  over  the  entire  lifetime  of  the 


Fig.  13.  Evolution  of  cumulative  bit  volume  under  MIP  with 
directional  antennas  for  D-MIP  and  RB-MIP;  T=  oo,  F  =  oo, 

a  =  4(P=l-,(pr,pR)  =  (0,0)). 


network  (in  this  case  until  no  pair  of  live  nodes  is  within 
communication  range),  as  a  function  of  t?m;n,  for  both  RB- 
MIP  and  D-MIP. 

Fig.  14a  shows  Bx,e  for  (pT,pR)  =  (0, 0)  and  (3  =  1. 
Consistent  with  the  results  presented  above,  D-BIP  provides 
better  performance  than  RB-BIP,  and  this  difference 
increases  as  f?m;n  decreases.  There  is  little  difference  in 
performance  for  6m\„  >  90.  However,  there  is  approximately 
an  order  of  magnitude  difference  for  0m;n  =  1  (the  smallest 
value  for  which  results  were  obtained). 

Figs.  14b  and  14c  show  the  impact  of  processing  power, 
for  the  cases  of  ( pT,pR )  =  (0.01,0.1)  and  (0.1,  1),  respec¬ 
tively.  The  most  obvious  impact  of  processing  power  is  the 
reduced  value  of  Bx,e-  Since  energy  is  now  expended  for 
signal  processing,  less  is  available  for  RF  transmission. 
Therefore,  the  overall  delivered  traffic  volume  is  reduced 
greatly  (note  that  the  vertical  scale  is  logarithmic  in  Fig.  14a 
and  linear  in  the  others).  Moreover,  the  advantage  of  using 
D-BIP  decreases  as  processing  power  increases,  again 
because  a  smaller  fraction  of  energy  is  available  for  RF 
transmission. 

8  Conclusions 

In  this  paper,  we  have  studied  the  impact  of  the  use  of 
directional  antennas  on  the  design  and  performance  of 
algorithms  for  energy-aware,  source-initiated,  session- 
based  broadcasting  and  multicasting  in  all-wireless  net¬ 
works.  Our  primary  focus  has  been  on  operation  in 
energy-limited  environments,  in  which  the  nodes  are 
subject  to  hard  constraints  on  available  energy.  By 
contrast,  in  applications  where  energy  efficiency  is  of 
primary  concern,  it  is  assumed  that  sufficient  energy  is 
available  to  support  all  communication  requirements,  but 
its  use  incurs  a  cost. 

The  main  contribution  of  this  paper  is  the  development 
of  four  tree-construction  algorithms  that  are  suitable  for  use 
with  directional  antennas,  two  for  broadcasting  and  two  for 
multicasting.  One  of  the  broadcast  algorithms,  Reduced- 
Beamwidth  BIP  (RB-BIP),  uses  the  trees  formed  by  BIP 
under  the  assumption  of  omnidirectional  antennas  and  then 
reduces  the  beamwidth  to  concentrate  the  RF  energy  in  the 
cone  where  it  is  needed.  The  other,  Directional-BIP  (D-BIP), 
exploits  the  directionality  of  the  antennas  throughout  the 
tree-construction  process.  The  corresponding  multicasting 
algorithms  are  RB-MIP  and  D-MIP.  The  Reduced-Beam- 
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Fig.  14.  Bit  volume  per  unit  energy  vs.  0min  for  D-MIP  and  RB-MIP  (/?  =  1).  (a)  (pT,pR)  =  (0,0),  (b)  (pT,pR)  =  (0.01,0.1),  and  (c)  (pT,pR)  =  (0.1, 1). 


width  approach  can  be  used  with  any  tree  construction 
algorithm  as  a  postalgorithm  step;  however,  the  Directional 
approach,  which  uses  the  incremental  power  philosophy 
we  developed  in  our  earlier  research  and  which  is 
characterized  by  higher  complexity,  is  integrated  into  each 
step  of  the  BIP  algorithm. 

We  have  shown  that  the  incorporation  of  residual  energy 
into  local  cost  metrics,  which  results  in  load  balancing  that 
spreads  the  burden  of  energy  use  among  more  of  the  nodes, 
has  a  considerable  impact  on  network  performance.  Most 
importantly,  we  have  shown  that  the  time  of  the  first  node's 
death  can  be  delayed  significantly,  thus  permitting  opera¬ 
tion  at  maximum  throughput  rates  much  longer  than  is 
possible  when  a  criterion  of  minimum-power  trees  is  used. 
Additionally,  the  overall  volume  of  data  that  is  delivered  is 
increased.  System  operation  is  highly  robust  with  respect  to 
the  residual-energy  parameter  (5;  values  between  0.5  and  2 
have  been  shown  to  work  well. 

The  tree  construction  algorithms  developed  in  this  paper 
for  directional  antennas  provide  significant  improvement  in 
terms  of  network  lifetime  and  total  delivered  traffic  volume, 
as  compared  to  their  omnidirectional  counterparts.  The 
improvement  is  greatest  for  small  values  of  9m jn  and  low  to 
moderate  values  of  signal-processing  power.  Moreover,  the 
“directional"  versions  provide  significantly  better  perfor¬ 
mance  than  the  "reduced-beam"  versions,  especially  for 
small  values  of  #m;n  and  small  values  of  processing  power. 
However,  they  do  so  at  the  cost  of  increased  complexity. 
The  relative  advantage  of  the  directional  versions  decreases 
as  the  propagation  constant  a  increases.  Thus,  the  simpler 


reduced-beam  versions  may  provide  acceptable  perfor¬ 
mance  in  many  practical  applications. 

Our  algorithms  are  heuristics,  and  no  claim  for  optim¬ 
ality  is  made.  Nevertheless,  by  illustrating  the  degree  of 
improvement  that  can  be  obtained  by  exploiting  the 
properties  of  directional  antennas  into  tree-construction 
algorithms,  we  have  demonstrated  the  benefit  of  the  use  of 
directional  antennas  in  energy-constrained  environments. 
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